Collective Modes of = 2 Quantum Hall Bilayers in Tilted Magnetic Field 
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We use time-dependent Hartree Fock approximation to study the collective mode spectra of = 2 
quantum Hall bilayers in tilted magnetic field allowing for charge imbalance as well as tunneling 
between the two layers. In a previous companion paper to this work, we studied the zero temperature 
global phase diagram of this system which was found to include a symmetric and ferromagnetic 
phases as well as a first order transition between two canted phases with spontaneously broken (7(1) 
symmetry. We further found that this first order transition line ends in a quantum critical point 
within the canted region. In the current work, we study the excitation spectra of all of these phases 
and pay particular attention to the behavior of the collective modes near the phase transitions. We 
find, most interestingly, that the first order transition between the two canted phases is signaled by 
a near softening of a magnetoroton minimum. Many of the collective mode features explored here 
should be accessible experimentally in light scattering experiments. 



I. INTRODUCTION 

Scattering experiments have provided extremely pow- 
erful and important probes of two dimensional electron 
systems^. A particularly nice application of light scat- 
tering was a recent set experiments'^''* on quantum Hall 
bilayers with equal densities in each layer. An apparent 
mode softening at total filling fraction ly = 2 was identi- 
fied with the existence of a Goldstone mode which fit well 
with prior predictions of a novel canted phase in these 
bilayer systems^. The experiments were conducted us- 
ing the tilted-field technique for sweeping across a wide 
range of Zeeman energies in situ. Interestingly, tilted 
magnetic fields also have a non-trivial effect on inter- 
layer tunneling. Furthermore, as noted first by Burkov 
and MacDonald^, in = 2 bilayer systems with a charge- 
imbalance between the layers, tilting the magnetic field 
can induce a new first order quantum phase transition 
embedded in the canted phase. As shown by the current 
authors in the preceding companion paper to this work^, 
the first order transition separates two phases with the 
same symmetry which are topologically connected in the 
phase diagram. Similar to a liquid-gas transition, the 
first order phase transition line terminates at a quan- 
tum critical point. In our preceding paper— we discussed 
the phases and phase transitions of = 2 in detail, ac- 
counting for both charge imbalance and in-plane mag- 
netic field. The purpose of the present work is to ex- 
amine the excitation spectra of these different phases in 
order to make connection with possible future experi- 
ments. Particular attention will be paid to the evolution 
of the collective-mode dispersions across the new first- 
order transition induced by the tilted magnetic field. 

Bilayer quantum Hall systems in general have been the 
focus of a great deal of recent study^. The already rich 
physics of quantum Hall effects is further enhanced in bi- 
layers by the added degree of freedom. The most studied 
of the bilayer quantum Hall state is certainly the v = I 
state^. At v = 1 the spin degrees of freedom are effec- 
tively frozen out, and all the interesting physics occurs 



in the isospin (layer degrees of freedom). In contrast, for 
i> — 2 systems, not only the layer but also the spin de- 
grees of freedom are important. Indeed, in several of the 
phases of v — 2, the spin and isospin degrees of freedom 
are actually entangled. 

In perpendicular magnetic field, v = 2 bilayers exhibit 
the many-body canted phase at finite tunneling-^. This 
phase is a spontaneously-broken U{1) symmetry phase, 
despite the finite tunneling. This is in marked contrast 
with the V — 1 bilayers, in which a U{1) symmetric 
phase is possible only in the absence of tunneling. Things 
change when a finite voltage bias is added: In this case, 
1^ — 2 bilayers can exhibit a many-body phase in the ab- 
sence of tunneling. This phase is somewhat akin to the 
many-body phase of = 1 bilayers, as was pointed out by 
MacDonald, Rajaraman, and Jungwirth^. These authors 
therefore mused that, in the presence of a finite in-plane 
magnetic field component, v = 2 bilayers may also un- 
dergo a commensurate incommensurate transition. 

Burkov and MacDonald^ explored this possibility. In- 
deed, they found that charge-unbalanced v = 2 bilayers 
can undergo a phase transition driven by the in-plane 
field component. However, the phase transition was be- 
tween two commensurate phases, instead of a between 
a commensurate and an incommensurate phase. One of 
the commensurate phases was akin to the commensu- 
rate phase of the v = 1 bilayers - the isospin compo- 
nent followed the magnetic field. The other commensu- 
rate phase however, was more peculiar: in this phase, 
both isospin and spin components were commensurate 
with the in-plane field! In our previous publication, we 
attempted to understand the physics behind this spin- 
commensuration. We explored the phase transition fur- 
ther, and found that it terminates at a critical point 
within the canted phase (See Fig.QJ. 

As mentioned above, one way of experimentally distin- 
guishing between the many phases of the v = 2 bilayers 
is to probe the collective excitations—. In this paper 
we therefore set out to explore this phase transition fur- 
ther by finding the collective modes. The many-body 
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phases (CI, and C2 as well as / which occurs in the 
absence of tunneling) of the u = 2 bilayers are charac- 
terized by spontaneously broken symmetries, which re- 
sult in the formation of low-energy Goldstone modes. In 
fact, the first theoretical and experimental evidence of 
the canted phase in charge-balanced v = 2 bilayers was 
obtained by observing a softening spin-density mode in 
time-dependent Hartree-Fock (TDHF) analysis^ and in 
inelastic light-scattering experiments2i4. 

More generally, the TDHF approximation allows one 
to predict the response of the system to any one of a 
number of possible experimental probesiSiii. A general 
perturbation can be introduced into the system by the 
addition of a small time-dependent term, 5H: 



to the total Hamiltonian (Eq. The operators p^cr.u 
are the density operators 



(1) 



3x 



-ik^X-ik^kyl'^ /2 J 



(where g is the Landau level degeneracy) whose time- 
evolution can be obtained using the Heisenberg equation 




Q, 



FIG. 1: Global phase diagram of charge-unbalanced v = 2 
bilayers in tilted magnetic field. The phase diagram is cal- 
culated for bare tunneling gap Ag^^ — 0.06e^/(eZ), Zeeman 
energy ~ 0.01e^/(eZ), and the distance between the lay- 
ers d = I. The axes are the amplitude of the in-plane field 



wavevector: Q|| 



{B±r^/d) 



and the external bias voltage A\ 



This choice of axes is particularly suitable, since current ex- 
perimental techniques allow to vary the bias voltage and the 
in-plane field m situ over a wide range of values. SC is the 
spin-singlet commensurate phase, CI is the simple, isospin- 
commensurate phase, C2 is the spin-isospin commensurate 
phase. Solid lines represent the second-order quantum phase 
transitions, dashed line — the novel first-order C1-C2 tran- 
sition. 



of motion (Eq. \27\i )'^'^: the external time-dependent field 
'l'po-,i/cr' (k, cj) in Eq. ^ must turn into its complex con- 



jugate when /id <^ va' and k, 



— k, —uj so that the 



Hamiltonian remains Hcrmitian. Different experimental 
probes will couple to different combinations of the den- 
sity matrix p^^o-.i/o-' (k). For example, a surface acoustic 
wave experiment might couple to the total charge density 
"^ua Pva,i^(j, whereas certain spin-polarized light scat- 
tering experiments might couple to the spin-flip-density 

If the perturbing external field ^^a,!^a' (k, uj) is small, 
one can assume that the response of the system to it is 
linear: 



,(k,w) $^^,i.^/(k,w). 



(3) 

The proportionality coefficient between the change in the 
density expectation value as a result of the perturbation, 
{Sp^'sys'iK^)) and the perturbing field, ^^a,,ya'{K^), 
is the density response function that can be obtained in 
the TDHF approximation. The presence of a pole at a 
particular frequency and wavevector indicates resonant 
response (i.e. the presence of a collective mode) . 

We therefore start the derivation of the collective-mode 
dispersion of the v ~ 2 bilayers by obtaining general ex- 
pressions for the density response function. First, in Sec- 
tion^ we review the unrestricted Hartree-Fock approxi- 
mation through which the ground states of the = 2 bi- 
layers were obtained in our previous companion paper». 
We then continue in Section IIIII to derive the TDHF 
approximationiSiiSii^ which is tailored to match the un- 
restricted Hartree-Fock of our prior study, and results in 
a general matrix equation for the density response func- 
tion. 

In Sections IIVI and IVl collective modes of the charge- 
unbalanced 1^ — 2 bilayers in perpendicular field are ob- 
tained (collective modes in charge-balanced 1^ — 2 bilay- 
ers were discussed in Refi^ and&). In Section Hvl the 
collective-mode dispersions (i.e. the poles of the density 
response function) of the charge-unbalanced ly = 2 bilay- 
ers with no interlayer tunneling are obtained in closed 
form. The symmetry properties that simplify the (com- 
plicated) general equations for the density response func- 
tion are discussed. In Section EI the collective- mode dis- 
persions of the charge-unbalanced v = 2 bilayers in the 
presence of a small amount of interlayer tunneling are ob- 
tained numerically and compared to the collective-mode 
dispersions of the system without tunneling (Figs.|3Q. 



Section IVII presents our main result: the collective- 
mode dispersions of the v = 2 bilayer systems in tilted 
magnetic field. A set of collective-mode dispersion curves 
calculated as the tilt-angle is increased and the system 
undergoes the C1-C2 transition is exhibited in Fig. 
FigureElshows the dispersion curves of the 1/^2 bilayers 
at the critical end-point. A dramatic softening of the 
Goldstone mode at this point is observed. 
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II. "UNRESTRICTED" HARTREE-FOCK 
APPROXIMATION: AN OVERVIEW 

We review here the "unrestricted" Hatree-Fock first 
discussed previously in Refs. and0 Our system con- 
sists of a disorderless zero temperature bilayer quantum 
Hall system with tunneling between the layers and both 
perpendicular and in-plane magnetic fields. Three terms 
of the Hamiltonian — Hz^ the Zeeman energy, Hy the 
bias voltage between layers, and Ht the tunneling — cou- 
ple to single electrons and comprise the non-interacting 
part of the Hamiltonian 



Ho = 



Hz + Hy + Hj^ 
X 



A 



SAS 



X 



where S and / are the spin and isospin operators 

fl ss' 
S (.LU 



(4) 

(5) 
(6) 



where a and r are sets of Pauli matrices. Here, the 
subscript X represents the momentum index of the elec- 
tron state in Landau gauge, and the subscript s takes 
on the values -f 1 and —1 corresponding to spin up and 
spin down whereas ^ and v take on the values -1-1 and 
— 1 corresponding to different layer index (up and down 
"isospin" ) . 

The Coulomb interactions between the electrons are 
taken into account by an additional term 

(72 

X c c 

'^Licri Xi+qyP /J2''"2 ^2 fi2'^2 X2+qyl^ /^icriXi 

where intra- and mterlayer Coulomb interactions are 



(7) 



VRR{q) 



2^ 

eq 



VRL{q) 



(8) 



respectively, d is the distance between the layers, and 
ri is the area of the sample. The total Hamiltonian is 
therefore simply 



H = Ho + Hi. 



(9) 



The Coulomb-interacting Hamiltonian in Eq. ^ is not 
tractable exactly, and we solve it using the Hartree-Fock 
approximation. In the usual manner, we assume that the 
many-body ground state, jG), is a Slater determinant of 



single-particle states and perform a functional minimiza- 
tion of the expectation value {G\H\G) with respect to 
these single-particle states. As was described in Refs. Q 
and '3, under the assumption of translational invariance 
in the y-direction, the trial ground state | G) can be writ- 
ten in the form 



|G) = 11/1^/2^10), 



where 



X 



fl(T 



-iQl^aX 



Ci_iaX ■ 



(10) 



(11) 



A ground state (Eq. with nonzero possesses 

spin-isospin-wave order, discussed at length in Ref. 0. 

As was mentioned in Ref. 's', the proposed ground state 
(Eq. and H10|l ) is not the most general Slater deter- 
minant (Hartree-Fock) state. However, our analysis of 
the collective modes around the ground states obtained 
by the minimization of {G\H\G) indicate the stability 
of these states against second-order transitions that can- 
not be described within the Hilbert space defined by our 
ansatz. We note that the possibility of phase transitions 
into a soliton-lattice state cannot be ruled out in this 
work. 

To obtain the approximate Hartree-Fock ground state 
for the v — 2 bilayer system, we minimize the expectation 
value of the Hamiltonian in Eq. 1^, ^{G\H\G) , with re- 
spect to the variational parameters z"^^, and Q^^o-. As was 
demonstrated in Ref. 0, the resulting set of minimization 
conditions can be arranged in the form of a Schrodinger 
equation: 



(12) 



where Z" — (z^f , zr^, z^f , z^i), and Af is a 4 x 4 matrix, 
which is just the mean-field single-particle Hartree-Fock 
Hamiltonian: 

M^s';t,s = -Az5f,„a^^^, ~ AvSss'T^,,- 

^ X 



i\2MiQ\\-Qi)x)ry,] 



lis s',m=l,2 

- F^.(-[Q//2 + Qsl2 (s - s')]q^) X 

X E i^^is'T^^- (13) 



n=l,2 



Here, a further simplification of the problem is made by 
making an assumption that 



Ilia — ^Qi + 2^3, 



(14) 



where a finite Qi indicates the presence of an isospin- 
wave order, while a finite Qs reflects the real spin- wave 
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order. The functions ff-(q) and Fp,y(q) used in Eq. H13() 
are defined as 



(27r)2 
dk 



iqAk( 



(15) 



27r 

1 1 - e^'''' 

(l^flfl(q) - Vrl{<1)) = -7 ^ , (16) 



el 2ql 



These functions arise from the Hatree and exchange parts 
of the interaction Hamiltonian (Eq. Q) treated in the 
Hartree-Fock approximation. 

The Schrodinger equation (Eq. (|12|l ) is solved 
iterativelji^. At each iteration, the two eigenstates corre- 
sponding to the lowest eigenvalues are filled (i.e. chosen 
to be the states 1 and 2). These lowest-energy eigenstates 
and Z'^ are then used to obtain the matrix M for the 
next iteration. The procedure is repeated until a self- 
consistent solution is achieved. This solution — a set of 
eigenspinors Z" sorted according to their eigenvalues — 
defines the lowest energy trial state among the Slater de- 
terminants defined by Eqs. l(Tn|l and ((TT|l subject to fixed 
values of the Q^s's- The eigenvalues e„ give the binding 
energy of a particle in the subband n, i.e. it is the en- 
ergy lost when the particle is taken out of the system. 
The sum of individual binding energies does not give the 
groundstate energy; the ground state energy is calculated 
from ^{G\H\G) (see Ref. 0). The minimization of the 
energy of the ground state over the Qfj.s's is done last^S,. 
Thus, we find the Hartree-Fock ground state in two steps: 
First, for fixed values of Qi and Qs, we minimize the ex- 
pectation value of the ground state energy with respect 
to (z^V)*- Then, we minimize the ground state energy 
with respect to Qj and Qg. Thus, we obtain the phase 
diagrams thouroghly discussed in Ref . (3, a representative 
example of which is given in Fig. ^ 



III. THE DENSITY RESPONSE FUNCTION 

We start our analysis of the response of the system to 
external perturbations by considering the possible exci- 
tations of electrons between the subbands. The presence 
of the layer and spin degrees of freedom of electrons in 
bilayer systems results in the splitting of each Landau 
level into four such subbands. When the filling fraction 
is ly ~ 2, the lowest two subbands are filled, i.e. the Fermi 
energy lies between the second and the third subbands. 
An elementary transition of a non-interacting i/ = 2 bi- 
layer system occurs when a particle is moved from one 
of the filled levels into one of the empty levels, resulting 
in a particle-hole pair. Four such transitions are possible 
in the v = 2 bilayers (provided the cyclotron energy is 
assumed to be much larger than all the other relevant 
energy scales): 1 ^ 3, 2 -> 4, 2 ^ 3, and 1 ^ 4. The 
energy of an unbound particle-hole pair is the difference 
between the energy gained by inserting the particle into 



an empty level, e^, where /3 = 3,4, and the energy lost 
by removing it from a filled level Eq., where a = 1,2: 
Ae — 6/3 — Eq. In a real system, the particles and the 
holes they leave behind interact. The interactions lower 
the energy of the particle-hole pairs and make it wavevec- 
tor dependent. The energies of the unbound particle- hole 
pairs show up as the poles of the HE density response 
function, which is, diagrammatically, the bare density re- 
sponse function dressed with self-energy corrections. The 
self-energy corrections represent the effect of the renor- 
malization of the single-particle levels by the interactions, 
accounted for in the HE approximation. In order to ac- 
count for the particle-hole interactions, the HE density 
response function is dressed with vertex corrections. In 
the charge-unbalanced i/ = 2 bilayers both the Hartree 
"bubbles" and the exchange "ladders" contribute to the 
collective-mode dispersions^. 

As was discussed in Ref. Isl we frame our problem so 
that the v — 2 bilayer ground state can have a very simple 
form, \G) — Hx fixAxl^)' ^^"^ basis of the creation- 
annihilation operators fnx (Eq. Ullfl ). It is therefore 
convenient to define generalized density operators 



^ X 



fa,xfp 



X+kyl^ 



(17) 



Note that the expectation values of the generalized den- 
sity operators are always diagonal in the ground state of 
the u = 2 bilayer systems. 



(18) 



Since in tilted magnetic fields, in the gauge of our choice, 
the operators contain the creation operators cj^^. with 
different position-dependent phase factors e*'^^"''^, the 
generalized density operators, Pap, are related to the 
physical density operators, p^a.ua' (Eq. |(5J), not only 
through a linear transformation but also by a shift of 
wavevector: 

p«^(k) ^ e-W-+Q.^-)^-«'^/^z^,(z,^,,)* X 

X Pp.a,iya'0^- {Qfia - Q^a')x). (19) 

Now, to perturb the system so as to determine its re- 
sponse, we rewrite the external perturbation Hamiltonian 
Eq. ^ in terms of the generalized density operators as 



5H 



where 



E 

7^ 



duj f d k ~ 



2tt / 2tt 



|.^5(k,c.)e-'"*p^,(k), (20) 



$^5(k,w)= 



X {QtL<y - Qua')x,Uj). (21) 

Thus, application of an external potential $^^,1^0-' , which 
couples to the physical density at wavevector k, can gen- 
erate perturbations $ that couple to the generalized den- 
sity at other wavevectors. For calculational simplicity. 
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and for simplicity of presenting our results, we will fo- 
cus on calculating the response of the system to $ which 
couples to the generalized density. From this result, one 
can simply determine the physical response of the sys- 
tem to an arbitrary perturbation (in terms of the phys- 
ical density). However, these wavevector shifts between 
the physical density and the generalized density must be 
kept in mind as they can be nontrivial, as we will see 
below rSec. rvrai and Eq. |73). 

To determine the response of the system to the time- 
dependent perturbation in Eq. (|20|l . we use standard lin- 
ear response theory (Kubo formula) , in which the result- 
ing change in the expectations of generalized density op- 
erator is assumed to be proportional to the perturbation: 

{5p^0{\^,u)) = Xo'/5V(k,^) %s{K^). (22) 

The proportionality coefficient is the retarded density re- 
sponse function 



xTp,s{K^)^^W e-*([p„Mk,i),p;5(k,0)]). (23) 
Jo 

We obtain the collective-mode dispersions of the v = 2 
bilayers by finding the poles of this response function. 

It is convenient to obtain the retarded density response 
function (Eq. (|23(l ^ from a corresponding imaginary-time 
density response function 



Xa/375(k,r) = -g{Tpal3{^^,T)p\s{\L,0)) 



(24) 



where pap = Pap — {Pap)- The imaginary-time density 
response function can be Matsubara-transformed to get 
Xct/37i5(k, *i^) (where iVLn are bosonic frequencies) that, 
in turn, can be transformed into the retarded density- 
response function x^^^^(k, cj) by a Wick rotation iVl 
LD -f iS. 

Following Cote and MacDonald (CM)ii, we proceed by 
calculating the Hartree-Fock density response function 

v° ( 



^(k, ifi) from its equation of motion: 



ld_ 
gdr 



Xa/37<5(k,T) 



(25) 



= 5(r)([p„^(k, 0), p;,(k, 0)]) + (T— p„0(k, T)p\,{k, 0)). 

The commutation relations of the generalized density op- 
erators paf3(k,0) are 

g[pa^(k), ^-(k)] = SpsPajiO) - SajPpsiO)- (26) 

The time-evolution of the density operator is determined, 
within the Hartree-Fock approximation, from the mean- 
field Hartree-Fock Hamiltonian. H^^: 



d 

— pap{\^,T) = [n"^ ,Pap{KT)] = 
OT 



(27) 



where the mean-field Hamiltonian is diagonal in the fa 
basis, and can simply be written as (see Sec. 0) 



(28) 



Using the Hartree-Fock equation of motion for the den- 
sity operators, and Matsubara-transforming the equation 
of motion for the density response function, we get 

in^xlp.siK^^n) ^ (29) 
= Spsipa-^iQ)) - Sa-iipspiQ)) + (e/3 -e„)x°/37^(k,^^^«)■ 
The single-particle density response function is therefore 

liln- ep + Ca 

Its poles are, as expected, at the single-electron excita- 
tion energies. 

To take into account the interactions between the 
single-particle excitations, following Refs. ITll andHol we 
introduce the vertex corrections. The vertex corrections 
in TDHF arc Hartree "bubbles" and exchange "ladders" , 
which have to be related to the Hartree-Fock self-energies 
through Ward identities^. Using the interaction con- 
stants 



a/375(k) — 2^^2 $Z ^iicri (^fiai) ^i'2CT2 (^»2CT2 ) 



ll ,l2.(Tl ,tT2 



il,i2,CTl,<T2 



X F^,^^{[h^ - {Q^,a, - Q»2<T2)]^ + kyV) (32) 

we dress the single-particle density response function to 
include the interactions between the single-particle exci- 
tations: 

Xa/375(k,il)„) = X°/375(k,iaO- 

- Xa/3ab(k, i^n)Fcabd(y)Xcd-iS{K i^n) (33) 
Xa/375(k, iri„) = Xa/375(k, ir2„)-|- 

+ Xa/3ab(k, i^n)Hbacd{k)Xcd-yS[k, i^n); (34) 

summation over repeated indices is implied, and the 
exchange interaction functions Fi^i^(\<i) are defined in 
Eq. 

To solve for the dressed density response function, 
Xap-y si)^, i^n), it is Convenient to cast Eqs. H29|) - (|34|l into 
matrix form. With the definitions 

Xia+p-SAy+S-a = XQ/37<5(k, ifi„), (35) 
^4q+/3-5, 47+5-5 = Sps{paj{0)) ~ Saj{p5p{0)) 

= Sa^ySpS E {Sa.n - 5p^n), (36) 
n=l,2 

M4^a+p-5Ay+5-5 = SaySps{ea — €p), (37) 

Hia+P-SAj+S-i — HpajsC^.), (38) 

F4a+p-5Al+S-5 — F^apsO^-), (39) 
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we have 



in„Xo = R-MXo, 
X = Xq — XqFX, 
X = X - XHX. 



(40) 
(41) 
(42) 



The density response function Xq/St* (k, if^n) is repre- 
sented by the 16 x 16 matrix X: 



X = [iVLn - M - R{H - F)]-^R. 



(43) 



The poles of the density response function arc the solu- 
tions to the secular equation 



det[iVLn -M - R{H - F)] = 0. 



(44) 



The 16 X 16 matrix equations can be reduced to 8 x 8 
by eliminating the forbidden single-particle excitations, 
such as 2 ^ 2, or 3 ^ 4. Even though the effective 
"energies" of these transitions are solutions to the secular 
equation (Eq. H44() ). it is easy to show that the weights 
of these modes are always and they do not show up 
in the density-response function matrix. The remaining 
8x8 matrix equation includes the interactions between 
the four transitions that create the particle-hole pairs, 
and four corresponding transitions that recombine the 
particles and the holes. 



IV. COLLECTIVE-MODE DISPERSIONS OF 
CHARGE-UNBALANCED v = 2 BILAYERS, 

AsAs = 

In the most general case, the solution to Eq. H44|) has to 
be found numerically. The situation simplifies consider- 
ably when either the bias voltage or the tunneling is zero. 
The former case has been studied in Refsi^ and&, who 
obtained the spin-density wave branches of the collective 
modes in perpendicular magnetic field (their analysis can 
be easily extended to the case of tilted fields) . The latter 
case is considered in this section: Using the symmetry 
properties of the v = 2 bilayers in the absence of tunnel- 
ing, we show that different inter-subband single-particle 
excitations are independent of each other m. v — 2 bi- 
layer systems in the absence of tunneling and the ver- 
tex corrections simply result in additional (g-dependent) 
renormalization of the excitations. 

In this section, we present the analytical calculation of 
the collective-mode dispersions in v = 2 bilayer systems 
in the absence of tunneling in perpendicular field. We 
explain the main features of the dispersion curves and, 
in the second part of this section, discuss the evolution 
of these features as the interlayer tunneling is turned on. 
The V — 2 bilayers in titled magnetic fields are considered 
in the next section. 



A. Parametrized Asas = ground state 

We start the calculation of the collective-mode disper- 
sions by finding the ground state of the system. As was 



discussed in Sec. the ground state of the v = 2 bilay- 
ers is obtained within the Hartree-Fock approximation 
by solving the Schrodinger-like equation (Eq. ((T^ ). In 
the absence of tunneling, the mean-field solutions Z" can 
be parametrized by two parameters, so that a transfor- 
mation matrix, S", that can be constructed of the four 
eigenspinors, Z", has the form: 



5= (z^z^z^z4) = 



/I 








sin 


e 


-i<i> Q 












COS! 



e*"^ COS 61 
-sine 




1/ 



(45) 

where Z" are defined after Eq. (|12|) . The two sub- 
bands with the lowest binding energies are filled. By 
construction, for positive bias voltage and Zeeman cou- 
pling, the lowest are the bands 1 and 2, so that the gen- 
eral form of the ground state is Y[ 



(sin 9 c\ 



Rix 



e'-'l' cos9 zl'\c]^^x)\^)- When cos 6' ~ 1, the ground state 
is the ferromagnetic state; when cos 6* = 0, it is the spin- 
singlet state; the intermediate values of cos 9 indicate 
that the system is in the many-body so-called /-state 
(See Ref. onlinecitelopatnikova:03a). It is easy to see 
that the Hamiltonian is invariant with respect to change 
of the phase cj) — this is the U{1) symmetry that results 
in the formation of a Goldstone mode in the /-phase. To 
simplify our calculations, we choose (j) — 0, so that the 
matrix S is now real and S = S~^. 

When the parameter 6 is such that the ground state 
energy is minimized, the S matrix diagonalizes the mean- 
field Hamiltonian matrix M (recall that the Schrodinger- 
like equation (Eq. H12I) ') is the result of a formal mini- 
mization of the Hartree-Fock ground state energy with 
respect to the parameters (z^a-)*)- We can therefore find 
9 by forcing the matrix A — S^^MS to be diagonal — 
i.e. equating the off-diagonal terms of the matrix A to 0. 
The resulting minimization condition can be written as 



Ko{9) sin 261 = 0. 
where the function Kq(9) is defined as 
Av-Az~ 2H^ 



(46) 



Ko{9) - - 



(//_ -F_)cos20. (47) 



Equation (|46|l is satisfied automatically in the ferromag- 
netic and spin-singlet phases, where cos^ — and 1, 
respectively, so that sin 26* — 0. In the /-phase, the equa- 
tion is solved by 



2 ^ 1 Ay - Az - 2//_ 

cos W = : ^ — 

2 4:{H^-F^) 



(48) 



where 



\2itV 



el2V 2 



(49) 
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and H- is defined in Eq. H16|) . Equation (|48() has a so- 
lution when its right-hand side takes a value between 
and 1. This region is the region of stability of the I- 
phase. (Note that, in the absence of anisotropy between 
the inter- and mtralayer interactions, Eq. (|48|l would have 
no solution, since = and F_ = in this case.) If 
the right-hand side of Eq. H48|l is negative, the ground 
state energy is minimized by cos^ 6 = and the system 
is in the ferromagnetic phase. If the right-hand side of 
Eq. 148|l is greater than 1 , then cos^ 9 = 1 and the sys- 
tem is in the spin-singlet phase. Intermediate values of 
cos^ 9 give the /-state. Note that the function Ko{9) = 
throughout the /-phase and takes on finite values in the 
ferromagnetic and spin-singlet phases. 

The resulting binding energies are the eigenvalues of 
the matrix A/, and can be read off the diagonal of the 
matrix A = S^^MS. 

ei = -Az + F^cos29 + Ko{d)-{F++F^) (50) 

€2 = F^ + Ko{9) COS 29 -{F+ + F^) (51) 

£3 = -F^ - Ko{9)cos29 (52) 

£4 = Az - F_cos29 - Koie). (53) 

Note that, in the many-body region, where Kq{9) = 0, 
the smallest single-particle gap depends only on the in- 
teraction constants, £3 — £2 = F^ — and is constant 
throughout the region. This is again a manifestation of 
the many-body nature of the /-phase. 



B. Symmetry properties of the Asas = ground 
states 

In Ref. we pointed out that the ground states real- 
ized in v = 2 bilayer systems in the absence of tunnel- 
ing and with positive Zeeman field and bias voltage are 
eigenstates of the operator /^ -I- with the eigenvalue 
g (where g is the Landau level degeneracy, i.e. -1-1 per 
flux quantum). The operator F -\- commutes with the 
bilayer Hamiltonian and therefore provides a good quan- 
tum number to classify the eigenstates of the Hamilto- 
nian. Thus, the ground state belongs to the class of states 
with the F -\- S'^-quantum number equal to g. So does 
the lowest-energy excited state, which is the result of the 
2^3 transition; that is to say, the transition 2-^3 does 
not change the F -\- 5'^-quantum number of the system, 
5{F + S"=) = 0. The excitations 1^3 and 2 ^ 4, on 
the other hand, lower the F- + S'^-quantum number by 
one, i.e. 5{I^- + S^) = —1, and the highest-energy ex- 
citation 1^4 lowers the F' -\- S'^-quantum number by 
5{F + S^) = -2. 

Modes characterized by the same 5{F + S^) can be 
mixed by the Coulomb interactions, unless they can be 
classified further by other quantum numbers. For the 
excitations 1^3 and 2 — > 4, simply the operator 



"AT 



Ci^i provides a good quantum number (it com- 
mutes with the zero-tunneling Hamiltonian): the exci- 
tation 2^4 possesses the same eigenvalue g as the 



ground state, while the excitation 1^3 changes it by 
Sn-fii = — 1. Therefore, all the four single-particle inter- 
subband modes are decoupled from each other as states 
with different conserved quantum numbers. Only the 
particle-hole interactions within the same mode therefore 
appear in the calculation of the collective modes, and ma- 
trix of the density response function thus separates into 
four 2x2 matrices. 



C. AsAS = collective-mode dispersions — general 

Reduced to include only one excitation mode a — *■ /3, 
and its counterpart (3 ^ a, the 2x2 density-response- 
function matrix obeys a matrix equation of the same gen- 
eral form as the full equation (Eq. (|43|l '). The matrices 
comprising the reduced Eq. H43|l are 



M=i '^'^''^ 







R = 



1 

-1 



The Hartree part of the vertex-correction matrix is 



H 



FlpaaP FlpafSa 





(Ho 


Ho \ 


)- 


[ Ho 


Ho j 



(54) 
(55) 

(56) 



where Ho = Hpaap and we used the symmetries of Hap-yS 
in Eq. H31|l with respect to the exchange of indices (given 
that, without loss of generality, the coefficients z'j'^ can 
be assumed to be real in the absence of tunneling). We 
use the same symmetries (Eq. (|32|l ) to obtain the general 
form of the exchange contribution to the vertex correc- 
tions: 



F = 



FaaPP FpaPa 

OLOi 



Fo Fi 
Fi Fq 



(57) 



where Fq = Faapp and Fi = Fapap- The solution to 
Eq. (|44|l gives the dispersion curve of the collective exci- 
tation a — > /3: 



Fo- Fi- [ep - £a)x 
-2Ho + Fo+Fi-{ep-e^) 



(58) 



where the £„ are given in Eqs. (jSOll - ltS^ . 

The resulting collective-mode dispersions are given in 
Figs.|2H3| 



D. Goldstone mode 



We start by considering the lowest-energy mode, which 
softens when the system enters the /-phase. Using the 
parametrization of the coefficients zf^ given by Eq. H45|l . 
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we get 



Ho = H3223 
Fq ~ -F2233 = 
Fi = -F2323 = 



= H- sin^ 29 
F. - cos^ 29 



sin^ 29. 

The resulting collective mode dispersion is 

co23 = VQi0,q) {e{e,q) + T{9,q)) 

where 



(59) 
(60) 
(61) 



(62) 



0(9, q) = 2X0(6*) COS 261 + F^{q) - F+iq) 
r{e,q)^2iH_{q)-F_{q))sm^29, 

and F^{q) = F„(g) - F„(0). (f±(q) - ^(i^flfl(q) ± 
-PitL(q)), where i^p,y(q) is defined in Eq. (fT5)l : i?-(q) is 
defined in Eq. ((TB|l 1 

The functions (q) are proportional to q^ at small val- 
ues of q, while approaches a finite constant. When 
the system is in the /-phase, where Ko{9) = and 
sin^ 29 ^ 0, the dispersion curve becomes gapless and 
<^23 oc l^l at small q. Indeed, this is the linearly dis- 
persing Goldstone mode that appears in the /-phase as 
a result of the spontaneously broken U{1) symmetry of 
the ground state. The Goldstone mode disperses linearly 
since the generator of symmetry does not commute with 
the Hamiltonian. The velocity of the Goldstone mode is 



VG 



= Zsin26'A/2(/f_ - F. 



d 



RL 



(63) 

where F^l = F+ - F_ is calculated in Eq. gSJl. The 
Goldstone-niode velocity is proportional to sin 29 — it is 
zero at the phase boundaries and the greatest near the 
middle of the /-phase. 

In the ferromagnetic and spin-singlet phases, the dis- 
persion curve of the mode is gapped and analytical 
around q = Q. In the ferromagnetic state the dispersion 



■^23 



2F_-Av + (F_(9)-F+(g)), 



(64) 



consistently with the fact that the ferromagnetic state is 
stabilized by the magnetic field and by the anisotropy of 
the Coulomb interaction in bilayer systems. As the bias 
voltage is increased, the gap is reduced, until it becomes 
zero when Ay = -I- 2F^. The transition to the /- 
state occurs at this point, consistently with Eq. (|48|) . The 
dispersion of the 2^3 mode in the spin-singlet phase is 



23 



A^ - 4//_ + 2F^ + (F^iq) - F+{q)). (65) 



In the absence of tunneling, the system is driven into 
the spin-singlet phase by the external bias voltage 
and against the renormalized interlayer charging energy 
2(2//_ -/"-). Again, consistently with Eq. the 
system undergoes a mode-softening phase transition be- 
tween the spin-singlet phase and the /-phase when Ay — 
Az-h2(2//_ -F_). 



The last interesting feature of the 2^3 mode we con- 
sider is the roton minimum that this mode develops in the 
/-phase (see Fig. ^ . The roton minimum appears deep 
in the /-phase and disappears close to the boundaries 
with the ferromagnetic and spin-singlet phases. It oc- 
curs at w 1 , a wavevector characteristic of interaction 
effectsi^. In the present case, the roton minimum indi- 
cates a tendency toward formation of an interlayer spin- 
density wave. (Formally, it is the non-trivial wavevector 
dependence of //_ (q) — F_ (q) in Eq. that causes the 
roton minimum to appear.) 



E. Spin-wave modes 

The dispersions of the modes 1 — + 3 and 2^4 have the 
form given by Eq. (|58|l . It is clear from the parametriza- 
tion of the zJJ^ (Eq. gSJ) that //„ = Fi = for the 
excitations 1^3 and 2 — > 4. The 2x2 interaction ma- 
trices H and F are therefore diagonal, and the dispersion 
curves have the simple form 



(66) 



Unlike the dispersion of the Goldstone mode 2^3, the 
dispersions of the 1^3 and 2^4 modes are analytical 
in all the phases of the = 2 bilaycrs in the absence of 
tunneling. The only relevant interaction constant, Fq, is 
the same for both modes, 1^3 and 2^4, 



Fq — Fi 



133 



F2: 



F+iq)-F^iq)cos29, (67) 



and so are the binding-energy differences, e^ — ei = ei—e2- 
The collective-mode dispersions are therefore degenerate, 
and 



^^^13 — ^24 



(68) 



= Az - F+{q) - F^{q) cos 29 - 2Ko [9) sin^ 9. 

The dispersion is always gapped; the gap is equal to 
the Zeeman splitting in the ferromagnetic and the /- 
phases. The system possesses a finite magnetization in 
these phases, so that the modes 1-^3 and 2 — > 4 corre- 
spond to spin-wave modes. In the spin-singlet phase the 
magnetization is zero, and the gap of 1 — > 3 and 2^4 
modes departs from the Zeeman splitting linearly with 
Ay: 

wfg = Lol^ = (69) 
= Ay - 4//_ -f 2F_ + (F_(g) - F+{q)). 

That the gap of cjfg is equal to the Zeeman splitting at 
the boundary of the spin-singlet and the /-phases is more 
clear if one compares this dispersion curve to w^g, given 
in Eq. H65|) . The difference between cjfg and wfg equals 
to the Zeeman splitting at any Ay and q. This is because 
the excited state produced by the excitation 2 — > 3 is a 
spin-triplet state with = -1-1. The excited states that 
result from 1 — > 3 and 2 — > 4 are superpositions of a spin- 
singlet state and a spin-triplet state with S"^ = 0. The 



excitation 1 — > 4, which we consider below, in the spin- 
singlet state results in a spin-triplet excited state with 
S'- = -1. 

The degeneracy of the modes 1^3 and 2 ^ 4 is a 
consequence of the up-down, left-right symmetry of the 
Hamiltonian that reverses the sequence of subbands: the 
transformation that exchanges the levels 1 and 4, and 2 
and 3 — ^ -Az, Av ^ -Ay, Asas ^ -^sas, 
and c^o- — > e^'^'^/^c-^ leaves the Hamiltonian invari- 
ant. This symmetry maps the transition 1 ^ 3 to 4 — > 2. 

F. Highest-energy mode 

For the highest-energy mode, 1^4, the interaction 
matrices also turn out to be diagonal, Hq = f i = in 
Eqs. l(5T|l and . The collective-mode dispersion there- 
fore has the same form as that for the 1^3 and 2^4 
modes (Eq. ^). The remaining exchange interaction 
constant, Fq, is 

Fo = F,,ii = F+{q)-F^{q), (70) 

and the dispersion relation, therefore, is 

a;i4 = Az + Ay-2(iJ_-i^_)- (71) 
- 2iJ_ cos 29 - F+ (q) + F. (q) . 

When Ay = 0, the mode 1 — > 4 is degenerate with the 
2 — > 3 mode (Eq. (|65|l ). The modes are split when a 
bias voltage is applied to the system, and the splitting 
grows as 2Ay until the system enters the /-phase (at the 
point where the 2^3 dispersion becomes gapless). The 
gap of the 1^4 mode starts decreasing as the system 
is brought deeper into the /-phase. At the boundary 
between the /-phase and the spin-singlet phase, the gap 
is at its lowest value of 2Az- In the spin-singlet phase, 
the dispersion is 

uf^ = Ay + Az - 4/r_ + 2F_ + (F^iq) - F+{q)), (72) 

and continues linearly with Ay, as it did in the ferromag- 
netic phase. The dispersion ujf^ is simply related and 
uji^, as the excitation resulting in the third {S^ — —1) of 
the triplet excited states. 

V. COLLECTIVE-MODE DISPERSIONS OF 
CHARGE-UNBALANCED = 2 BILAYERS, A / 0, 

The Hamiltonian of the charge-unbalanced i/ = 2 bi- 
layer systems in the presence of interlayer tunneling does 
not commute with the /^ -I- operator. The symmetry 
considerations that we used to obtain an analytical solu- 
tion for the density-response function of the v = 2 bilay- 
ers in the absence of tunneling cannot be used to simplify 
Eq. (|44|l when tunneling is present. We therefore use nu- 
merical techniques to calculate the dispersion relations 
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FIG. 2: The collective-mode dispersions of the charge- 
unbalanced 1/ — 2 bilayers in perpendicular field - ferromag- 
netic phase. The dispersions in the left column are given for 
a system with Ag^g = 0.0(e^/eZ); in the right column - for 
a system with Ag^g = Q.Q2{e^ / el); the Zeeman energy in all 
panels is Az = 0.01(e'^/eZ); the bias voltage is given in the 
upper right corner of each panel in units of {e^ /el). Some 
collective-mode dispersions are degenerate. The low-energy, 
low-wavevector region is shown in the insets. 



of the charge-unbalanced v ~ 2 bilayers with finite inter- 
layer tunneling. Much insight into the numerical results 
can be gained by comparing the collective-mode disper- 
sions in the systems with tunneling to the analytical re- 
sults for the systems without tunneling. The comparison 
is presented in Figs. [3 El and 01 for the ferromagnetic, 
many-body, and spin-singlet phases, respectively. The 
sets of plots in each figure are arranged in two columns: 
In the left column, we plot the dispersion curves of a 
system with no tunneling. The dispersions of a system 
with tunneling are given in the right column. In all the 
plots, the Zeeman energy is set at A^ = 0.01(e^/eZ); the 
external bias voltage is given in the upper right corner of 
each panel. 
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FIG. 3: The collective-mode dispersions of the charge- 
unbalanced V = 2 bilayers in perpendicular field - many-body 
phase. The dispersions in the left column are given for a 
system with Ag^^^ = Q.Q{e^ / el); in the right column - for a 
system with Ag^g = 0.02(e^/ei); the Zeeman energy in all 
panels is Az = 0.01(e'^/£/); the bias voltage is given in the 
upper right corner of each panel in units of {e^ /el). Some 
collective-mode dispersions are degenerate. The low-energy, 
low-wavevector region is shown in the insets. 



A. Ferromagnetic phase 

In Figure 13 we present the dispersion curves obtained 
in the ferromagnetic phase. While the ferromagnetic 



ground state does not change as the bias voltage is in- 
creased, the collective-mode dispersions demonstrate the 
evolution of the inter-subband energetics that eventually 
leads to a phase transition. The top two panels of Fig. |21 
show the dispersion curves in the absence of bias volt- 
age. In the absence of tunneling, we can identify the 
lower curve as the degenerate 1 — > 3 and 2^4 exci- 
tations, that have a gap equal to the Zeeman energy as 
g — > 0. The Zeeman gap is hard to discern in the figure, 
since scale of the Zeeman energy is very small in compar- 
ison with the energy scales of the other excitations. The 
upper curve represents the excitations 1 — > 4 and 2-^-3 
that are degenerate in the absence of tunneling and bias 
voltage. When a small amount of tunneling is present, 
the degeneracy of the dispersion curves is lifted, and the 
curves are split by an energy of order 0{/!^'gj^g/ F-) — a 
minute energy difference, not visible on the scale of the 
figure. 

As the bias voltage is increased, as was discussed in the 
previous section, in the absence if tunneling, the split- 
ting between the dispersions of the 1^4 and 2 — > 3 
modes increases linearly with the voltage. Thus, at 
Ay = 0.3(e^/£Z), the splitting of the upper branches is 
very large and apparent. The Zeeman branch, 1 — > 3 and 
2 ^ 4, in the systems without tunneling, remains degen- 
erate for any bias voltage. Note that, in the absence of 
tunneling, the dispersion curve of the 2 — s- 3 mode crosses 
the Zeeman branch clearly without interacting with it. 
The situation changes when tunneling is present: the (ap- 
proximately) 2 — > 3 mode develops an anti-crossing with 
one of the modes of the Zeeman branch. The interactions 
between different excitations are allowed in the presence 
of tunneling by the broken symmetry of the Hamilto- 
nian: when the tunneling term is present, the operator 

-I- does not commute with the Hamiltonian; the 
eigenvalues of the operator are, therefore, no longer good 
quantum numbers of the excited states. Nevertheless, 
one mode always stays independent of other excitations: 
One can see that, in the states with finite magnetiza- 
tion. Figs. 121 and 01 the spin- wave mode, identified by 
the Larmor minimum at the Zeeman energy, is always 
decoupled from the other modes. This is a consequence 
of the up-down, left-right symmetry of the Hamiltonian 
mentioned in the Sec. IIV El and preserved in the pres- 
ence of tunneling. This symmetry maps the excitation 
1 ^ 3 to 4 — > 2. The mapping results in a special form 
of the density-response matrix that always separates out 
one Zeeman mode. 



B. Many-body canted and /-phases 

As the voltage is increased further, the gap of the 2 — > 3 
mode decreases until it disappears. The softening of the 
mode signals the onset of a many-body phase — the I- 
phase in the system without tunneling and the canted 
phase in the presence of tunneling. The critical voltage 
is higher in the system without tunneling, since, when 
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FIG. 4: The collective-mode dispersions of the charge- 
unbalanced V = 2 bilayers in perpendicular field - spin-singlet 
phase. The dispersions in the left column are given for a sys- 
tem with Ag^5 — 0.0{e^ /el); in the right column - for a 
system with A^^g = 0.02{e^ /el); the Zeeman energy in all 
panels is = 0.01(e^/eZ); the bias voltage is given in the 
upper right corner of each panel in units of {e^ /el). Some 
collective-mode dispersions are degenerate. The low-energy, 
low-wavevector region is shown in the insets. 



the velocity soon reaches nearly a constant, while in the 
/-phase it continues increasing until approximately the 
middle of the /-phase. As the system gets closer to the 
transition to the spin-singlet phase, the Goldstone-mode 
velocity goes to in the reverse fashion. 

Another effect of the broken-symmetry of the Hamilto- 
nian in the presence of tunneling is the further widening 
of the gap that develops in the mode that splits off the 
Zeeman branch as a consequence of the mixing with the 
other modes. This gap is much larger than /S.sas and 
is therefore due to the interactions. Finite Ag^g in this 
case mainly serves to break the symmetry of the Hamil- 
tonian. This, "third" , mode develops a non-analyticity 
a,t q — Q and a ring of shallow roton minima, degenerate 
for all directions of q. 

The roton minimum that, as we showed in the pre- 
vious subsection, characterizes the /-phase is inherited 
by the low-tunneling canted phase. It is however less 
deep in the canted phase, and gradually disappears as 
tunneling is increased until it becomes more important 
than the bias voltage (i.e. when the bias voltage does 
not result in a significant charge imbalance within the 
many-body phase). Another feature to observe is the 
lowering of the energy of the highest-energy mode. Its en- 
ergy scale changes dramatically as one sweeps across the 
many-body phase by increasing the bias voltage. Around 
the boundary between the ferromagnetic phase and the 
many-body phase, the gap of the highest-energy mode is 
two orders of magnitude larger than (it is rather of 
©(Ay)), but at the boundary of the many-body phases 
and the spin-singlet phase, the gap decreases to 2Az. 



the tunneling is present, it works in concert with the 
bias voltage to stabilize the spin-singlet phase and desta- 
bilize the ferromagnetic phase. The details of the col- 
lective mode dispersions as the gap of the 2^3 mode 
approaches zero around Ay « 0.6(e^/e/) are given in 
the insets of the bottom two panels of Fig. |21 In both 
insets, one can clearly see the Zeeman branch with the 
gap at the Zeeman energy, A^ = 0.01(e^/eZ). In the ab- 
sence of tunneling, Zeeman branch is degenerate; when 
tunneling is finite, the degeneracy is lifted by the inter- 
action with the 2 — > 3 branch near q = Q. The inter- 
action of the 2 — > 3 excitation with the superposition 
of the 1 — > 3 and 2^4 excitations result in an ear- 
lier onset of the many-body phase and a mixed mode 
with a large gap around q = Q. When the bias voltage 
is increased a little above Ay w 0.6(e^/£Z), the system 
undergoes a phase transition into the many-body phase. 
The collective-mode dispersions for the v — 2 bilayers in 
the many-body phase are given in Fig. |31 Right after the 
transition the velocity of the Goldstone mode is very low, 
but it rapidly increases as the system is taken deeper into 
the many-body phase by an increasing bias voltage. Near 
the transition, the velocity of the mode increases faster 
in the canted phase, than in the /-phase. In the canted 
phase, however, because of the inter-mode interactions. 



C. Spin-singlet phase 

When the bias voltage is around Ay w 1.4(e^/e/), the 
system undergoes a phase transition from the many-body 
to the spin-singlet phase. At this point the Goldstone 
branch develops a gap. The Goldstone mode becomes 
the lowest of the three spin-triplet excitations above the 
spin-singlet ground state. As can be seen in the insets, 
the three modes are separated by A^. In the absence 
of tunneling the S*^ = spin triplet is degenerate with 
the only allowed spin-singlet excitation. When tunnel- 
ing is present, the energy of the spin-singlet excitation 
is affected by the interactions; it slowly approaches the 
energy of the S*^ = spin triplet as the increasing bias 
voltage turns the interlayer-phase coherent spin-singlet 
state at Asa5 into a, v ^ 2 spin-unpolarized mono- 
layer state. 



VI. COLLECTIVE-MODE DISPERSIONS OF 
CHARGE-UNBALANCED v = 2 BILAYERS IN 
TILTED FIELD 

As we discussed in Sec. above and in Ref. '3, tilting 
the magnetic field away from the normal to the plane of 
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FIG. 5: The collective-mode dispersions of the charge- 
unbalanced V = 2 bilayers in tilted field - evolution across 
the C1-C2 phase transition. The dispersions are given for 
a system with A| = 0.01(e^/e;), A?;^^ = 0.06(e^/d), and 
Ay = 0.8(e^/e/). The strength of the in-plane component of 
the tilted magnetic field is given by the wavevector Q\\, which, 
for every set of dispersion curves, is represented by a number 
in the units of Xjl. When Q\\l < 0.96, the system is in the CI 
phase; when Q|[/ > 0.96, it is in the C2 phase. In the bottom 
left panel, the dotted line is for comparison of the dispersion 
curves in the C2 phase to those in the /-phase. 



the bilayer system leads to new phases and, in charge- 
unbalanced systems, phase transitions. In the charge- 
balanced ly = 2 bilayers, the tilted fields do not change 
the topology of the phase diagram. The interlayer phase 
coherent phases - spin-singlet and canted - become com- 
mensurate, and the ferromagnetic state is not affected by 



the in-plane field. In commensurate interlayer phase co- 
herent phases, the interlayer exchange interactions are ef- 
fectively destabilized by the in-plane field, and the phase- 
space volume of the canted phase decreases as the mag- 
netic field is tilted. The weakening of interlayer exchange 
interactions renormalizes collective mode dispersions of 
the charge-balanced v = 2 bilayers but does not result in 
interesting new features. 

The presence of a finite in-plane component of the mag- 
netic field produces more interesting effects in charge- 
unbalanced systems. When tunneling is strong enough, 
the in-plane field induces a new phase transition between 
the simple commensurate phase, CI, stable at low in- 
planc fields, and the spin-isospin commensurate phase, 
C2, more favorable at higher in-plane fields. As is shown 
in Fig. Qand discussed in depth in Ref. 8,, the phases CI 
and C2 are connected to each other, and the first-order 
transition between them terminates at a critical point. 
To further study this first-order transition, we obtain a 
series of collective-mode dispersions, calculated for the 
CI and C2 states as the in-plane field is increased. We 
also obtain the collective-mode dispersions at the critical 
point terminating the C1-C2 transition. 



A. Collective-mode dispersions across the CI— C2 
phase transition 

The evolution of the dispersion curves within the 
canted phases as the magnetic field is tilted is given in 
Fig. |31 We choose a system with = 0.01(e^/eZ) and 
Ag^g = 0.06(e^/£/), and hold the external bias voltage 
at Ay = 0.8{e-^/el), so that the system is approximately 
in the middle of the canted phase (c.f. Figs.P). For each 
probed point on the phase diagram in Fig. ^ we plot 
in Fig. 13 the cross-sections of the collective-mode disper- 
sions in two directions - perpendicularly to the in-plane 
field (in the x direction in our calculations), and in the 
direction parallel to it (the y direction); these plots are 
given side by side. The measure of the magnitude of the 
in-plane component of the magnetic field, the wavevector 
is given as the number in the left panel. 

The top two panels show the collective-mode disper- 
sions in perpendicular magnetic field. In perpendicular 
field, the dispersions are the same in all directions, so the 
dispersion curves in the left and the right panels coincide. 
One can see the features discussed for the canted phase 
in the charge-unbalanced = 2 bilayers fSec. llVjl : the 
linearly dispersing Goldstone mode that has a roton min- 
imum around ql ~ I, characteristic of charge- unbalanced 
systems; the spin-wave mode that decouples from the 
other three modes and has a gap equal to the Zeeman 
energy (the resolution of the figure does not allow to see 
the Zeeman splitting because of the relatively small Zee- 
man energy); the large interaction-induced gap at q = 
of the third mode. The highest-energy mode is not visible 
in this figure. 

When the magnetic field is tilted, the collective modes 
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start changing: they become asymmetric with respect to 
q± —q±- The velocity of the Goldstone-mode in the 
negative -direction becomes greater than that in the 
positive gj^-direction. The roton minima also become 
asymmetric — they develop a lowest point in the nega- 
tive qx-direction. This behavior is reminiscent of the be- 
havior of the collective-mode dispersions of a (charged) 
superfluid under the influence of an external electromag- 
netic gauge field, A. In a superfluid, much like in the 
v — 2 bilayers in the canted phase, a U(l) symmetry 
is spontaneously broken. The symmetry breaking results 
in the formation of a linearly dispersing Goldstone mode. 
When an external field is applied to a charged superfluid, 
the superfluid order parameter acquires a twist, e^^^'^, 
and the Goldstone mode acquires an anisotropy: 

= |kJ«2-i( — )2|AP- . k, (73) 
V 2 mc mc 

where vq is the initial velocity of the Goldstone mode. 
While the in-plane field does not couple to the symmetry- 
breaking order parameters in = 2 bilayers in the same 
way it does in a superfluid, it does result in winding phase 
factors e*'3s^, e*'^^^, and e*(<3s±Q/)x_ Because of the 
gauge symmetry of our system, an equivalent picture can 
be drawn up, in which fictitious gauge fields proportional 
to Qsj Qi, and Qs ± Qi couple to the corresponding 
(uniform) order parameters. In fact, this is exactly what 
we have done, when we chose to work in the basis of 
the creation-annihilation operators (see Sec. IIIII and 
Eq. H17|ll. in terms of which the ground state is uniform. 

Because the in-plane field generates different phase fac- 
tors (e'<3s^^ f,^QIX^ g»(gs±Q/)^) for different or- 
der parameters, the v = 2 bilayer system is somewhat 
more complicated than a model superfluid. It is clear 
in Fig. |31 that different collective-mode branches (and 
different parts of the branches) do not respond to the 
presence of an in-plane field in the same way. Thus, the 
Goldstone mode "tilts" to the right in Fig. |S1 while the 
roton minima "tilt" to the left. An effective theory of 
coupled superfluids that would explain the behavior of 
the collective modes, as well as the C1-C2 superfluid- 
superfluid phase transition is a potential direction of fu- 
ture research. 

In addition to the anisotropics, the application of the 
in-plane field results in an increased Zeeman mode gap, 

= A^y^l + Q'^^l'^/d?. It is also apparent, especially 

at higher in-plane field, that the minimum of the Zee- 
man mode shifts from gj^ = to gj^ = Qs- This ef- 
fect is also a consequence of the our choice of the gauge. 
As was mentioned in Sec. IIIII the density response func- 
tion Xa;a7i5(k;T) that we calculate is gauge-dependent, 
because it is the response of the system to the excitation 
with the density operator /9a^(k). This density opera- 
tor is related to the physical density operators not only 
through a linear transformation, but also through a shift 
of the wavevector (Eq. H19(l 'l. A real excitation there- 
fore will pick out signals from different dispersion curve 



branches (according to allowed symmetries) at different 
wavevectors. Thus, for example, the physical response 
function for a real spin-flip excitation of wavevector k is 
a superposition of the Xap-ts{(\) t) at wavevector h + Qsx 

XMiMT;MiMT(k,^) = -.9(7'PMiMT(k,T)p|,x^|(k,0)) = (74) 

The dispersion curves in the direction parallel to the in- 
plane field stay symmetric with respect to q\\ ^ —q\\. As 
the tilt-angle is increased, they change mainly because 
the minimum of the Zeeman and other branches shift 
away from the q± = plane plotted in the right column 
of Fig. 

As the magnetic field is tilted further, the anisotropy 
of the Goldstone-mode velocity becomes greater, and 
the energy of its roton minimum around q±l ~ —1 de- 
creases. Near the Cl~C2 phase transition, = 0.96, 
the roton minimum becomes lower than the Zeeman en- 
ergy. However, before it reaches zero and the system be- 
comes unstable, the phase transition to the C2 phase oc- 
curs, marked by an abrupt change in the collective-mode 
dispersions. The entire picture is effectively shifted by 
Qsc2 ~ Qsci (s'^e Fig. El in the positive q± direction. A 
Goldstone mode appears in place of the roton minimum, 
and a roton-minimum replaces the Goldstone mode. The 
minimum of the Zeeman mode jumps from Qsci to Qsc2 j 
as expected from Eq. l(7^ and Fig. El 

As the in-plane field is increased further, the roton 
minimum of the Goldstone branch becomes less deep - it 
approaches the Zeeman energy, Az, from below. When 
the in-plane component of the magnetic field becomes 
large, Qn ~ 2, the C2 phase becomes very close to an 
/-phase. This is reflected in the collective mode disper- 
sions. In the bottom panels of Fig.jSJ at we can see 
the symmetry of the Goldstone branch reappearing (the 
dotted line in Fig.|31is given as a guide to the eye). One 
of the Zeeman branches of the /-phase forms an anti- 
crossing with the Goldstone branch at q± ~ 1. When 
Q|| « 2 the Zeeman modes and the Goldstone branch 
become nearly decoupled. The Zeeman branches in our 
calculation have a minimum at Qs ~ Q\\ in high mag- 
netic fields. This is again the gauge effect we described 
above. Zeeman modes shifted by —Qs, as they would be 
in terms of physical densities, would restore the /-phase- 
like appearance of the collective-mode dispersions deep 
in the C2 phase. 

B. Collective- mode dispersions at the critical point 

In the last part of this section, we consider the col- 
lective mode dispersions at the critical point. Again, 
we use a bilayer system with = 0.01 (e^/e/) and 
Ag^g = 0.06(e^/e/). The critical point for this sample 
occurs at Ay = 0.5328 and Q|| = 0.9362; Qs « 0.474 for 
these parameters. The spin- wave wavevector Q5 is hard 
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FIG. 6: The collective-mode dispersions of the charge-unbalanced v = 2 bilayers in tilted field at a critical point terminating 
the C1-C2 transition. The dispersions are given for a system with = 0.01(e^/e/) and Ag^5 — 0.06(e^/eZ). The critical 
point for this system occurs at Av = 0.5328 and Q\\ — 0.9362; Qs ~ 0.474 for these parameters. The top plot zooms in on the 
long-wavelength region, to demonstrate the dependence of the Goldstone-mode dispersions in the positive q± direction. 



to define precisely at the critical point, since the energy 
profile as a function of Qs is very flat: {Qs — QscY- 
The "flatness" of the energy as a function of the spin- 
wave wavevector Qs implies the existence of very soft 
spin-wave fluctuations. Indeed, as shown in Fig. [SJ the 
velocity of the Goldstone mode in the positive q± di- 
rection becomes (and the next-order in q± emerges: 
uj cx (7^). 

VII. CONCLUSIONS 

In summary, we have obtained the collective mode dis- 
persions of the charge-unbalanced v = 2 bilayers in tilted 
magnetic fields using the time-dependent Hartree-Fock 
approximation (Figs. [^HHI)- The collective modes pos- 
sess a number of characteristics that can be observed in 
light scattering experiments — such as softening modes 
and roton minima. Thus, the novel phase C1-C2 phase 
transition, discussed at length in RefH, is signaled by a 
near-softening of a roton minimum. 

The collective-mode dispersions of the 1^ — 2 bilayers 
in tilted fields exhibit behavior suggestive of a system 
of coupled superfluids under the influence of an external 
gauge field. Thus, when the magnetic field is tilted, i.e. a 
finite in-plane magnetic field is added to the system, the 



modes become Doppler-shifted (Eq. H73|l ). The Doppler 
shift is varies for different modes, reflecting on the fact 
that the in-plane magnetic field couples to the order pa- 
rameters of the v = 2 bilayer systems in different ways 
(see Sec. lvra . 

An interesting future direction of research therefore 
would be the construction of an effective model with two 
order parameters that spontaneously break U{1) symme- 
try. An external gauge field can couple to the order pa- 
rameters differently, so that when the breakdown of one 
superfluid occurs, the other superfluid is stable. Such a 
superfiuid transition would be an interesting model for 
the C1-C2 transition in the ly = 2 bilayers. 
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Note that the density operators p^^.^^/ (k, a;) are closely 
related to the (total) physical density n(k) as n(k) = 
^g-fc I II p,iCT,^CT(k, 0), where g is the Landau level de- 
generacy. 



